Plant Molecular Biology

Iron (Fe) homeostasis is fundamental to plant growth and development and is centrally regulated by a well-conserved network of basic helix-loop-helix (bHLH) transcription factors. In Arabidopsis thaliana, bHLH121 regulates the expression of FER-like iron deficiency-induced transcription factor (FIT) and clade Ib bHLH38/39/100/101 genes, among others. However, its functional conservation remains unexplored, and some contradictory results have been reported regarding the effects of bHLH121 overexpression. Here, we investigated the evolutionary conservation of bHLH121 through complementary tests in the Arabidopsis bhlh121 mutant created by CRISPR/Cas9 gene editing. We also investigated the effects of ectopically expressing full-length, N-terminal, and C-terminal fragments of bHLH121 on plant responses to Fe deficiency in Arabidopsis. Transgenic plants overexpressing the full-length bHLH121 driven by the cauliflower mosaic virus 35S promoter (bHLH121 ox) resulted in shorter roots and more severe leaf chlorosis compared to wild-type plants under Fe limiting conditions. Moreover, overexpression of the N-terminal fragment (N121 ox) rendered plants highly sensitive to Fe deficiency, whereas overexpression of the C terminal fragment (C121 ox) did not. Consistent with the phenotypic observations, N121 ox plants accumulated lower Fe content, and the ferric chelate reductase activity was significantly reduced. In N121 ox plants, the induction of key Fe regulatory genes, especially FIT, was significantly inhibited under Fe deficiency. Confocal imaging showed that both N121-GFP and C121-GFP fusion proteins could localize to the nucleus. Electrophoretic mobility shift assay showed that both full-length bHLH121 and N121 could effectively bind the probes derived from the promoter of FIT. The present work reveals a conserved function of bHLH121 in Fe homeostasis across several dicot species and demonstrates that the N121 may possess a dominant nature in Fe deficiency responses.

Amaranthus tricolor is a climate-resilient orphan crop whose betalain pigments contribute to nutritional quality and stress tolerance. However, breeding efforts have traditionally relied on field yield trials due to the limited availability of genomic information and genetic resources. To address this gap, we first used the Illumina and PacBio sequencing platforms to assemble a novel reference genome of A. tricolor, using the Japanese cultivar Bayam. By combining this reference genome with a core collection developed in our previous study, we aimed to identify genomic regions associated with pigmentation-related traits across winter and summer seasonal environments. Betalains are antioxidant phytonutrients that contribute to abiotic stress tolerance, making the regulation of their biosynthesis under heat stress particularly relevant for developing climate-resilient crops such as A. tricolor. As a result, ATR1.0CH03_5131103, which is associated with ATR1.0ch03g000565, was identified under the winter seasonal environment. Meanwhile, transcriptomic profiles in this study revealed that ATR1.0ch03g000565 was differentially expressed under both control and heat stress conditions, suggesting that it may be associated with the regulation of CYP76AD1 expression under optimal conditions. Under the elevated-temperature treatment, transcript levels of CYP76AD1 and 5GT decreased in Am087, coinciding with lower betacyanin-dominated absorbance indices. cDOPA-like homologs were upregulated; however, pigmentation-related absorbance remained reduced. Further molecular experiments, such as CRISPR-Cas9 knockouts or overexpression studies, are necessary to validate the regulatory mechanisms of CYP76AD1 and 5GT under heat stress.

BackgroundN6-Methyladenosine (m6A) regulates RNA fate in flowering plants, but the repertoire and expression patterns of m6A-related genes remain incompletely characterized in liverworts. Here, we characterized these gene families in the liverwort Marchantia polymorpha within a comparative framework spanning diverse plant lineages. ResultsA comparative survey of 119 plant genomes using sequence similarity and conserved-domain criteria identified 18 candidate m6A-related loci in M. polymorpha, comprising seven writer-related genes, five AlkB homolog genes, and six YTH- or CPSF30-related genes. Across the surveyed genomes, canonical writer components generally had low copy numbers, whereas AlkB homolog and reader-related families showed greater copy-number variation. Phylogenetic and domain analyses supported the family assignments of the M. polymorpha candidates. Eight of the 18 candidate loci showed detectable collinearity with Marchantia quadrata. Putative promoter analysis identified motifs annotated as hormone-, light-, and stress-responsive, with hormone-responsive motifs accounting for the largest share of predicted sites. Public transcriptome analysis revealed variation in candidate gene expression across developmental stages and tissues. Quantitative reverse-transcription PCR showed higher transcript abundance for all six YTH- or CPSF30-related genes in thalli than in gemmae or gemma cups. Cold, heat, osmotic, salt, and dark treatments elicited gene-specific expression responses. Darkness increased transcript abundance of several candidates, including MpMTA, MpALKBH1A, and MpECT2. A predicted protein-association network suggested links between subsets of the candidates and RNA-processing proteins. ConclusionsThese findings define a candidate m6A-related gene repertoire in M. polymorpha and identify tissue- and stress-responsive genes for targeted functional studies of RNA methylation in this liverwort.

High temperatures during grain filling stage are becoming increasingly frequent, compromising both grain and eating quality and thereby driving demand for heat-resilient cultivars. Such conditions are known to reduce the expression of granule-bound starch synthase I (GBSSI) and starch branching enzyme IIb (BEIIb), which are involved in starch biosynthesis, resulting in a decrease in amylose content and an increase in long-chain amylopectin. Thus, the present study introduced functional mutation in starch synthase IIa (SSIIa) that increases the proportion of short amylopectin chains to genetically compensate for the high-temperature-induced increase in amylopectin long chain. Rice lines carrying the ss2a mutation were grown at two locations with cooler (Akita) and warmer (Okayama) temperatures. Their grain traits, starch structure, and eating quality were compared. The ss2a mutant lines showed an increased proportion of short amylopectin chains as well as an increased apparent amylose content. Furthermore, these alterations in starch structure varied with the grain-filling temperature of the cultivation sites, ultimately affected eating quality. These results suggest that enriching short amylopectin chain via the ss2a mutation can counteract the increase in long amylopectin chain caused by high temperatures during grain filling, thereby maintaining a desirable starch structure and eating quality.

Understanding the genetic architecture of ecophysiological traits related to carbon gain and water loss is crucial to develop cultivars tolerant to abiotic stresses and that maintain productivity in the context of climate change. This can be achieved both through experiments in the greenhouse on potted plants under controlled conditions and in the field. However, the stability of genetic determinism of ecophysiological traits across environmental conditions has rarely been evaluated. We studied a grapevine diversity panel grown both in a greenhouse and in the field, to investigate the impact of growing conditions on leaf ecophysiological traits related to leaf morphology, fluorescence, gas exchange and water use efficiency. The correlation structure between traits was conserved within each environment even though, for most traits, cultivars ranking was not conserved between both environments. However, a set of cultivars, sharing the same functioning pattern in terms of gas exchange and fluorescence in both environments, was identified. Single- and multi-environment Genome-Wide Association Studies (GWAS), based on a dense genotyping matrix, were used to dissect the genetic architecture of leaf traits. The Quantitative Trait Loci (QTLs) identified through single-environment analyses differed substantially from those detected using multi-environment ones. Multi-environment approaches highlighted shared QTLs across environments for gas exchange and fluorescence traits. Candidate genes associated with carbon- and water-related functions were further identified in these genomic regions. Despite the weak correlation between greenhouse and field trait values, combining both environments enabled the identification of robust genetic regions associated with ecophysiological traits of interest.

Cucurbit yellow vine disease (CYVD) is a phloem disease caused by a group of Serratia ureilytica bacteria within the S. marcescens complex. The CYVD pathogen is distinct from other phloem pathogens in two fundamental ways: i) unlike other bacterial phloem pathogens, the CYVD pathogen readily grows on standard laboratory media; and ii) whereas most bacterial phloem diseases are transmitted solely by hemipteran insects, CYVD can be transmitted by both hemipteran (squash bug) and non-hemipteran (cucumber beetle) vectors. The phloem of cucurbits is also distinct, with cucurbits having a unique secondary phloem system, the extrafascicular phloem (EFP), for which a role in CYVD has not been investigated. The distinctive features of this pathosystem prompted a re-evaluation of pathogen localization in planta using a green fluorescent protein-tagged CYVD strain during infection of squash (Cucurbita pepo). In epifluorescence and confocal microscopy images of transverse and longitudinal sections of stem and petiole tissues, the fluorescence patterns of the CYVD pathogen were consistent with pathogen colonization of all four known EFP tissues, namely the peripheral, entocyclic, ectocyclic, and commissural sieve tubes. The CYVD pathogen may enter the EFP during herbivory by pathogen-infected vectors that cause mechanical damage, as the EFP exhibits much slower blockage following injury than the conventional fascicular (bundle) phloem. Moreover, although EFP has been proposed to have an antimicrobial function, EFP exudate collected from the stem and petioles of C. pepo supported rapid pathogen growth. These results provide evidence for a unique tissue tropism of the newly emerged plant pathogen S. ureilytica.

Drought stress severely limits tree growth and survival, a problem intensifying under climate change due to the rising frequency of extreme weather events. Identifying genetic variants that modulate drought stress responses and enhance resilience is therefore of great interest for developing climate-adapted trees. While the MADS-box transcription factor AGL16 has been implicated in drought responses and stomatal regulation in herbaceous species, the function of its homolog in woody perennials remains largely unexplored. Here, we investigate the role of PtaAGL16 in the poplar hybrid Populus tremula x P. alba. We generated CRISPR/Cas9 knockout lines of PtaAGL16 in the poplar clone INRA 717-1B4. Plants were subjected to a severe drought stress regime followed by a re-watering recovery phase. Morphological and physiological traits, including growth, stomatal conductance, shoot vitality, and stomatal density, together with stress-associated metabolites and drought-responsive gene expression, were analyzed. All agl16 knockout lines exhibited significantly reduced drought-induced damage compared to the wild-type. While the wild-type suffered extensive defoliation and persistent stem necrosis, the knockout lines retained higher shoot vitality and resumed growth after re-watering, including the formation of new leaves. Consequently, vitality in the knockout lines reached up to 86%, compared to only 34% in the wild-type. This improved resilience was accompanied by reduced persistent oxidative membrane damage, as indicated by lower MDA levels after recovery. These findings indicate that the disruption of PtaAGL16 is associated with enhanced drought resilience in a woody perennial species and highlight its potential as a genetic target for improving climate adaptation in fast-growing forest trees.

Boron toxicity constrains grapevine (Vitis) production in irrigated arid and semi-arid regions, yet most available rootstocks are sensitive to excessive boron. Although underutilized Vitis species may harbor sources of adaptation, the genetic basis of natural variation in boron exclusion remains poorly understood, limiting the development of boron-tolerant rootstocks. We screened 267 wild Vitis accessions representing 16 species for their ability to limit leaf boron accumulation under excess boron conditions. Leaf boron concentration varied 190-fold among accessions and was highly heritable (entry-mean H^2 = 0.71). Its correlation with stomatal conductance was modest (r = 0.34), indicating that transpiration rate explained only a small proportion of the variation in leaf boron accumulation. Leaf boron was correlated with leaf chloride (r = 0.36), but their largely non-overlapping association profiles supported distinct genetic control. Genome-wide association analysis identified a suggestive association on chromosome 11. Independent quantitative trait locus mapping in a full-sib F1 population (n = 158) derived from Longii-9018 (V. acerifolia) x GRN3 (commercial rootstock) identified qBEx11.1 in the same syntenic region on chromosome 11, which explained 19.5% of the phenotypic variance. Positional candidate genes near the association peak included a pectinesterase and a MATE efflux transporter, suggesting roles for cell wall modification and membrane transport. Leaf spectral reflectance was associated with leaf boron concentration in both populations, particularly at visible and red-edge wavelengths, indicating potential for high-throughput phenotyping. Together, these findings identify qBEx11.1 as a priority locus for boron exclusion and support integrated genetic and spectral approaches for breeding boron-tolerant grapevine rootstocks.

Background: Plant defensins are small cysteine-rich antimicrobial peptides. In Arabidopsis thaliana, the inducible Class I gene AtPDF1.2a is widely used as a transcriptional marker of jasmonate/ethylene immunity, whereas the constitutively expressed Class II members have received less functional attention. AtPDF2.2 was compared with AtPDF1.2a during infection of Arabidopsis by the vascular fungal pathogen Verticillium longisporum, the necrotrophic fungus Sclerotinia sclerotiorum and the bacterium Pseudomonas syringae pv. tomato DC3000. Results: Profiling of the AtPDF family showed that AtPDF2.2, together with the Class II genes AtPDF2.3 and AtPDF2.5, was significantly reduced during V. longisporum infection, whereas AtPDF1.2a was induced. A promoter-reporter line showed a corresponding decrease in AtPDF2.2 promoter activity during V. longisporum infection, strongest in infected roots. KD-pdf2.2 plants contained more V. longisporum DNA, whereas OE-PDF2.2 plants contained less. The same reciprocal disease pattern was observed during S. sclerotiorum and Pst DC3000 infection. Histochemical Prom::AtPDF2.2::GUS reporter signal decreased during S. sclerotiorum infection but increased following Pst DC3000 infection. Basal AtPDF2.2 expression was lower in JA and ET-signaling mutants and higher in the SA-deficient sid2 background, while the analyzed promoter lacked the GCC-box that mediates inducible AtPDF1.2a regulation. AtPDF1.2a overexpression reduced symptoms across the three pathosystems, but loss of AtPDF1.2a increased disease only during V. longisporum infection under the conditions tested and did not alter Pst DC3000 DNA abundance. AtPDF2.2 overexpression was not associated with reduced vegetative growth under the conditions tested, whereas AtPDF1.2a overexpression reduced rosette size and leaf number. Conclusions: AtPDF2.2 is a constitutively expressed Class II defensin whose expression level is associated with disease outcome across two fungal pathogens and a bacterium. Its promoter shows tissue-dependent regulation during V. longisporum infection, while histochemical reporter signal decreases during S. sclerotiorum infection and increases during Pst DC3000 infection. This regulatory behavior differs from that of the inducible reference gene AtPDF1.2a. The reciprocal disease phenotypes of the AtPDF2.2 transgenic backgrounds support a protective association that may involve direct antimicrobial activity, modulation of host responses or both. The absence of reduced vegetative growth in OE-PDF2.2 under the conditions tested supports further evaluation of PDF2-type genes in Brassica.

Cannabis sativa L. cultivation in tropical regions or under indoor conditions requires supplemental lighting to maintain plant yields. Night-break (NB) photoperiods offer a cost-effective alternative to conventional day extension, although genotypic responses have not been clearly described. This study evaluated the effects of continuous and pulsed NB treatments on the growth, phenology, physiology, productivity, and cannabinoid accumulation of three Cannabis sativa genotypes representing two chemotypes. Plants were exposed to a commercial 18-h photoperiod (T1), a continuous 1-h NB (T2), or four 15-min NB pulses (T3) with equal total duration. Continuous NB produced responses comparable to the commercial photoperiod. In contrast, pulsed NB accelerated phenological progression by two to three weeks and reduced final plant height by 15-28%, depending on genotype, without impairing photosynthetic performance or decreasing floral dry biomass. Harvest index increased under both NB treatments, indicating enhanced biomass allocation to reproductive structures, while cannabinoid accumulation remained primarily genotype-dependent with only minor treatment-specific modulation. These findings demonstrate that NB efficacy depends not only on the duration of nocturnal light but also on its temporal distribution, highlighting temporal organization of nocturnal light as a key regulator of photoperiodic responses and a promising strategy for optimizing cannabis cultivation.

Gene regulatory networks (GRNs) establish the spatio-temporal gene expression patterns necessary for reproducible development of organisms. MicroRNAs have emerged as important regulators of gene expression controlling either target gene levels, spatial patterns or timing. These molecular functions may vary depending on how they are integrated into GRNs. Here we used the Arabidopsis leaf as a model to study the genetic and molecular interactions between the CUP-SHAPED COTYLEDON 2 (CUC2) transcription factor and its negative regulator, miR164A. Quantitative morphometric analyses showed that miR164A acts from early stages during leaf development regulating CUC2 expression dynamic rather than the spatial distribution of CUC2. At the molecular level, CUC2 directly activates MIR164A expression, forming a single negative feedback loop. Further, using a ratiometric approach, we identified a cis-regulatory element (CRE) sufficient to drive CUC2-dependent MIR164A expression at the leaf sinus. Ordinary differential equation (ODE) modelling suggests that CUC2-CRE may control the timely accumulation of CUC2. Accordingly, quantification of CUC2 protein accumulation dynamics using a spatio-temporal mean cartography strategy in CRISPR cis-regulatory variants at the MIR164A locus revealed that CUC2-dependent MIRA164A expression is required to maintain both the precise timing of CUC2 expression and CUC2 levels. Altogether, our work provides experimental evidence for a specific role of the CUC2-MIR164A module in shaping leaf development, highlighting the contribution of network topology to GRN functions.

Abstract Ribosome profiling is widely used to map translation across the genome, yet current sORF discovery pipelines are constrained by their reliance on annotated protein coding regions and by treating biological replicates primarily as tools for quantification rather than as strict filters for reproducibility. This results in two major issues, the systematic exclusion of pseudogene derived loci from analysis and the generation of extremely large candidate sets in which true translation events are difficult to distinguish from noise. To overcome these limitations, a five-stage framework was developed that expands the search space to include pseudogene loci while enforcing dual biological replicate concordance as a hard binary filter. When applied to four human cell lines with matched replicates and benchmarked against existing tools, the framework produces dramatically smaller and more reliable candidate sets, reducing thousands to millions of predicted sORFs down to a small number of high confidence concordant loci. Across HeLa cell cycle conditions, a limited set of unique concordant sORFs was identified spanning multiple transcript classes including upstream ORFs, retained introns, and nonsense mediated decay transcripts, along with a pseudogene derived candidate. Notably, the pseudogene encoded sORF RPL26P19 was detected consistently in HeLa cells but missed or inconsistently identified by existing tools, while a separate pseudogene sORF from UBE3AP2 was observed specifically in normal primary cell types and absent in HeLa. Both candidates were later found to be supported by entries in OpenProt, even though this was not known during the analysis, which gives added confidence that the results are genuine and not biased by prior expectations. Overall, these findings suggest that pseudogenes are not just inactive remnants in the genome but can act as real, functional sources of translated sORFs in a cell type specific manner. At the same time, it becomes clear that most current tools are not designed to detect these signals effectively. This highlights the need for approaches that both broaden the regions being searched and apply stricter reproducibility checks, so that truly meaningful translation events can be identified more reliably.

Considering the increasing water and labor constraints in rice production, direct-seeded rice is becoming increasingly important. However, direct-seeded rice is vulnerable to water stress and inefficient nutrient use. Nitrogen use efficiency in rice is generally low, and improving nitrogen acquisition and assimilation is important for sustainable productivity. In this study, 48 recombinant inbred lines were evaluated under ammonium and nitrate nutrition in irrigated and rainfed environments across two cropping seasons. Agronomic, physiological, biochemical and molecular responses were assessed. Significant effects of genotype, nitrogen form, environments and their interactions were observed for several traits. Rainfed conditions reduced grain and biological yield compared with irrigated conditions. Ammonium produced a higher mean grain yield than nitrate by approximately 26.2 g m-2 under irrigated conditions and 154.2 g m-2 under rainfed conditions, although the responses varied among genotypes. Biological yield and SPAD values were generally higher with ammonium, whereas total chlorophyll showed season and environment-dependent responses. Glutamine synthetase activity was highest under ammonium treatment, while glutamate synthase activity was higher under nitrate treatment. Grain yield showed strong positive correlations with biological yield, but correlations between grain yield and GS or GOGAT activity were not significant. Expression of OsGln1.2, OsGln2 and OsGlt1 varied according to genotype, tissue and nitrogen treatment. Overall, the results demonstrate that responses to nitrogen forms are genotype- and environment-dependent, and that ammonium may provide a yield advantage under rainfed conditions. Integrating field performance with biochemical and molecular traits may help identify rice genotypes suited to water-limited systems.

DNA methylation dynamics regulate diverse plant stress responses, but the mechanisms governing arsenic (As) tolerance through epigenetic pathways remain largely unknown. Here we report that Arabidopsis mutants lacking the DNA demethylase, REPRESSOR OF SILENCING (ROS1), show substantially enhanced As resistance. Integrating transcriptome and methylome data, we found that PIN-FORMED 2 (PIN2), encoding an auxin efflux carrier that also transports As out of cells, represents a key ROS1 target. In ros1, we observed increased methylation in the last exon of PIN2, a region functioning as the promoter for an antisense long non-coding RNA (lncRNA). Hypermethylation of the promoter suppresses the antisense lncRNA, which follows reduction of sense PIN2 transcription. As a result, PIN2 expression increases in ros1 mutants during As exposure, enhancing As tolerance. Analysis of ros1 pin2 double mutants confirmed that PIN2 acts downstream of ROS1 in As tolerance. Our study demonstrates an epigenetic pathway in which ROS1 regulates As tolerance by controlling PIN2 expression through antisense lncRNA-mediated transcriptional interference, offering new epigenetic strategies using plants in heavy metal contaminated environments.

Chlorella sorokiniana is an industrially important microalga, but its genetic manipulation remains limited by poor transformation efficiency. Although several barriers have been proposed, the role of intracellular nuclease activity has received little attention. Here, we compared the degradation of exogenous double-stranded DNA in whole-cell lysates from transformation-recalcitrant (Chlorella sorokiniana and Tetradesmus obliquus) and transformation-proficient (Phaeodactylum tricornutum, Microchloropsis gaditana, and Nannochloropsis oceanica) microalgae. C. sorokiniana rapidly degraded both linear and plasmid DNA, while T. obliquus exhibited intermediate activity. In contrast, no detectable degradation was observed in the transformation-proficient species. DNA degradation in C. sorokiniana occurred within minutes, and nuclease activity in both recalcitrant species was inhibited by EDTA and Zn2+, consistent with a metal ion-dependent endonuclease. These findings suggest that intracellular nuclease activity may represent an underexplored barrier to exogenous DNA stability during transformation and provide a potential target for improving genetic transformation of recalcitrant microalgae.

Understanding how grafts form has fascinated scientists for over a century, and numerous transcription regulators involved in wound healing, callus formation, and grafting have been identified in model plants. However, the genetic architecture underlying grafting traits remains poorly characterized. To address this gap, we aimed to identify quantitative trait loci (QTLs) associated with grafting success rates and grafting-related traits (including callus formation, mechanical strength of the graft union, and post-grafting growth) in a grapevine rootstock mapping population, Vitis berlandieri x V. rupestris. Although no QTLs were identified for grafting success rate, QTLs were identified for component traits revealing a genetic architecture that would have remained hidden when considering grafting success as a single composite phenotype. By integrating QTL and transcriptome data, we identified regulators potentially mediating grafting-related traits. We then showed that one transcription factor that regulates callus formation and grafting in Arabidopsis potentially has a similar role in grapevine by silencing with exogenous double-stranded RNAs and that exogenous application of methyl jasmonate increases callus development at the graft interface. Candidate transcriptional regulators located within QTLs and evidence from gene silencing and hormone treatments suggest the conservation of regulatory mechanisms underlying wound-induced callus formation across plant species.

The apoplast is a primary location for interactions between plants and invasive pathogens and is the site where many pathogen-secreted effectors are perceived by cell surface immune receptors to activate defence responses. However, our understanding of apoplastic interactions remains limited because protein interactions are difficult to investigate in this harsh extracellular environment. TurboID-based proximity labelling (PL) has emerged as a powerful approach for studying protein interactions in plants, though its application has to date been restricted to intracellular proteins. Here, we designed and validated a TurboID-based PL strategy for investigating protein interactions in the leaf apoplast using the well-characterised interaction between the Phytophthora infestans elicitor INF1 and the receptor-like protein (RLP) REL in Nicotiana benthamiana. Transient expression of SP-INF1-TurboID (INF1-T) induced a cell death (CD) response comparable to that triggered by native INF1, demonstrating that fusion of the TurboID tag did not impair INF1 recognition by REL. Apoplastic localisation of both INF1-T and the control construct SP-eGFP-TurboID (eGFP-T) was confirmed, validating their suitability for PL experiments. Efficient TurboID-mediated biotinylation was achieved in the apoplast using co-infiltration of biotin, ATP and magnesium acetate. Streptavidin-HRP immunoblotting revealed distinct biotinylation profiles for INF1-T and eGFP-T. Furthermore, co-immunoprecipitation demonstrated specific biotinylation of REL by INF1-T, but not by eGFP-T, in wild-type, bak1 and sobir1/sobir1-like N. benthamiana. These findings demonstrate that TurboID-based PL is functional in the apoplast and provides a proof-of-concept for investigating elicitor-receptor interactions in this compartment.

Canonical ribonucleases such as Argonaute (AGO) participate in RNA-based gene regulatory processes. They are important for small RNA (sRNA) association and the formation of RNA-induced RNA silencing complexes. Potato (Solanum tuberosum) has 14 AGO-encoding genes, including the Solanum-specific StAGO15, which is activated upon infection by filamentous pathogens. The resistant phenotype conferred by ectopic expression of StAGO15 was abolished in plants in which a 21-nt sequence in the PAZ-encoding region of StAGO15 was targeted using an artificial miRNA. We then analyzed the sRNAs associated with StAGO15 in Phytophthora infestans leaf infection. Small RNAs were coimmunoprecipitated, subjected to deep sequencing, and analyzed against the P. infestans and potato genomes, in which cleavage sites were predicted by applying the R package, smartPARE. The most notable difference between pathogen- and water-inoculated StAGO15 samples was observed in the transposon (TE) category. Here, two categories of nonautonomous TE dominated: the miniature inverted-repeat transposable elements (MITEs) and short interspersed nuclear elements. The StAGO15-associated sRNAs bind and target complementary endogenous potato RNAs. One example is siRNA122, which acts on a cleavage site in the resistance gene PGSC0003DMG400015681 (StRES1). siRNA122 is derived from a MITE; upon deletion of this source sequence, resistance is restored. Our findings provide new insights into the plant immune surveillance network and demonstrate that deleting a 21 nt RNA sequence from a MITE transposon released a defense response, against two important filamentous pathogens.

Chlorophyll fluorescence imaging is a powerful tool used to characterize abiotic and biotic stress response of photosynthesis in plants. However, extracting metrics such as Fv/Fm (maximum quantum yield of photosystem II) from chlorophyll fluorescence images remains a substantial bottleneck in large measurement campaigns, often requiring significant manual effort. To resolve this issue, we developed fvfmPy, a user-friendly program which enables rapid analysis of chlorophyll fluorescence images of arrays of leaf discs. fvfmPy uses algorithmic image segmentation to identify regions of interest and extract fluorescence metrics from tens or hundreds of leaf samples per image with little to no manual input. fvfmPy is developed in Python and includes a graphical user interface enabling intuitive parameter adjustment and real-time quality control checking. We additionally developed fvfmR, an R wrapper package for fvfmPy. Fv/Fm values extracted with fvfmPy were indistinguishable from those obtained using manual methods (r2 > 0.99) but required 94% less time to complete (13 s vs. 200 s per 100 measurements). fvfmPy and fvfmR were built to directly integrate with existing fluorescence software packages and can greatly increase throughput of fluorescence measurements to help address key challenges in agriculture and ecology.

Premise:: Conifers comprise one-sixth of global forests and landscape genomic studies are urgently needed to understand their biology, yet relatively few conifer species have reference genomes. Targeted exon-capture is a practical alternative to whole-genome sequencing, but selecting loci without a reference genome remains a barrier. We developed a multi-tissue reference transcriptome for pi non pine (Pinus edulis) to ask how multi-tissue sampling improves probe design, and what is lost under a conventional single-tissue approach. Methods: We assembled a reference transcriptome from four diploid tissues (needle, phloem, root, vegetative bud) and one haploid tissue (megagametophyte), using the latter as a pseudo-guide. We designed a target-capture probe-set by prioritizing highly expressed P. edulis transcripts and transcripts matching a conifer exon database. Results: We recovered 76% of conserved plant genes with few duplicated transcripts. Filtering yielded 15,664 transcripts totaling 31.4 Mb, ~0.1% of the genome. Tissue-specific expression was substantial: 27.4% of transcripts occurred exclusively in non-needle tissues. Functional annotation confirmed complementary contributions, including seed-desiccation and germination genes from megagametophytes, and defense- and stress-associated genes from roots and phloem. Discussion: We built a transcriptome and efficient probe set using a roadmap suitable for other non-model plants with large genomes; keys to success were prioritizing complementary tissues and using expression-based filtering.